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ABSTRACT 

We study the effect of contamination by interlopers in kinematic samples of galaxy 
clusters. We demonstrate that without the proper removal of interlopers the inferred 
parameters of the mass distribution in the cluster are strongly biased towards higher 
mass and lower concentration. The interlopers are removed using two procedures pre- 
viously shown to work most efficiently on simulated data. One is based on using the 
virial mass estimator and calculating the maximum velocity available to cluster mem- 
bers and the other relies on the ratio of the virial and projected mass estimators. We 
illustrate the performance of the methods in detail using the example of A576, a cluster 
with a strong uniform background contamination, and compare the case of A576 to 15 
other clusters with different degree of contamination. We model the velocity dispersion 
and kurtosis profiles obtained for the cleaned data samples of these clusters solving 
the Jeans equations to estimate the mass, concentration and anisotropy parameter. 
We present the mass-concentration relation for the total sample of 22 clusters. 

Key words: galaxies: clusters: general - galaxies: clusters: individual: A576 - galax- 
ies: kinematics and dynamics - cosmology: dark matter 



1 INTRODUCTION 



The problem of contamination of kinematic samples of 
galaxies in clusters by foreground and background galax- 
ies has been recognized since such data became available. It 
arises because of the fact that only the projected positions 
and velocities of galaxies are measured in redshift surveys. 
Due to the lack of knowledge about the motion perpendic- 
ular to the line of sight it is difficult to judge a priori which 
of the galaxies found close to the cluster in projected space 
are actually bound to it and good tracers of the underlying 
potential. Including unbound galaxies, or interlopers, in the 
samples used for dynamical modelling may lead to signifi- 
cant bias in the estimated parameters of mass distribution 
in clusters. 

Many procedures have been proposed in the literature 
to deal with this problem. All these procedures aim at clean- 
ing the galaxy sample from non-members before attempting 
the proper dynamical analysis of the cluster. In the case of 
single clusters such a direct approach is the only one possi- 
ble because of small samples. Only in the case of kinematic 
samples combined from the data for many clusters one can 
attempt to take the presence of interlopers into account sta- 
tistically (van der Marel et al. 2000; Mahdavi & Geller 2004). 

The traditional direct approach by Yahil & Vidal (1977) 
relies on calculation of the velocity dispersion of the galaxy 
sample a and iterative removal of outliers with velocities 



larger than 3(t. In rich enough samples one can take into 
account the dependence of a on the projected distance from 
the cluster centre R and perform the rejection procedure in 
bins with different ct or fit a simple solution of the Jeans 
equation to the measured cr[R) profile and reject galaxies 
outside the 3a"(-R) lines (Lokas et al. 2006a). Perea, del Olmo 
& Moles (1990) discussed another method relying on itera- 
tive removal of galaxies whose absence in the sample causes 
the biggest change in the mass estimator. One can also com- 
bine the information on the position and velocity of a galaxy 
to infer the maximum velocity available to cluster members 
as proposed by den Hartog & Katgert (1996). 

The methods listed above (or variations of them) were 
recently tested in detail by Wojtak et al. (2006) using the 
results of cosmological Af-body simulations. Working with 
mock data samples created from the simulated dark mat- 
ter haloes they verified what fraction of unbound particles 
is actually removed from the sample. They found that the 
methods of den Hartog & Katgert (1996) and Perea et al. 
(1990) are the most efficient. 

Lokas et al. (2006a) have recently modelled the kine- 
matics of six nearby regular galaxy clusters with very few 
interlopers which could be reliably removed with the sim- 
ple 3(t(-R) scheme. Modelling both the velocity dispersion 
and kurtosis profiles allowed them to estimate not only the 
mass and concentration parameter of the cluster, but also 
the anisotropy of galactic orbits. On the opposite end, in 
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terms of the difficulty in estimating the cluster properties, 
lies Af689. As discussed by Lokas et al. (2006b) the cluster 
velocity distribution is so complex that no reliable estimates 
of its mass and concentration can be obtained from the kine- 
matics. 

In this paper we study the kinematics of clusters which 
could be classified as an intermediate case between the two 
mentioned above. We focus on significantly contaminated 
clusters, for which more sophisticated methods of interloper 
removal have to be applied, but still reliable estimates of the 
parameters can be obtained. We illustrate the performance 
of these methods in detail using the example of a rich nearby 
{z = 0.04) galaxy cluster Abell 576 (section 2). The cluster 
has a regular image in X-rays which suggests that it is in dy- 
namical equilibrium. However, contrary to similarly regular 
clusters like those studied by Lokas et al. (2006a) it is rather 
strongly contaminated by foreground and background galax- 
ies which are hard to separate from the main body of the 
cluster. We show how the inferred parameters of the cluster 
mass distribution would be biased if no interloper removal 
scheme was applied. Then we discuss the performance of dif- 
ferent interloper removal procedures. Finally we model the 
velocity dispersion and kurtosis profiles of the cleaned sam- 
ples to estimate reliably the virial mass, concentration and 
anisotropy parameter of the cluster. 

In order to place the results for A576 in the right context 
we study the kinematics of 15 more nearby clusters showing 
different degree of contamination (section 3). We demon- 
strate that after careful removal of interlopers such clusters 
can be modelled using Jeans formalism and present the re- 
sults of fitting their velocity dispersion and kurtosis profiles 
in terms of estimated mass, concentration and anisotropy 
parameter (section 4). This allows us to extend the mass- 
concentration relation obtained from the kinematic data by 
Lokas et al. (2006a) to the total of 22 clusters. 



2 MODELLING OF A576 

2.1 Contamination of the data by interlopers 

The positions and redshifts of galaxies in the vicinity of 
A576 were obtained from the NASA/IPAC Extragalactic 
Database (NED). We have searched the database for galax- 
ies within projected distance of 3 Mpc from the cluster cen- 
tre and with cz velocities differing from the cluster mean 
z = 0.0389 by less than ±4000 km s'^ The redshifts of 
galaxies thus selected come mainly from Mohr et al. (1996) 
and Rines et al. (2000). The cluster possesses two central 
galaxies. We choose one of these galaxies with redshift closer 
to the cluster mean as the centre of the cluster, with respect 
to which all distances are measured. 

The presence of this central galaxy pair may suggest 
that the cluster underwent a merger. This hypothesis is also 
supported by X-ray data: although the large scale X-ray sur- 
face brightness distribution measured by Einstein (Mohr et 
al. 1996) and ASCA satellites (Horner et al. 2000) is rather 
regular, some departures from equilibrium in the cluster core 
have been detected with Chandra (Kempner & David 2004; 
Dupke & Bregman 2005). However, since the X-ray map is 
not strongly perturbed and the centre of the gas distribu- 
tion coincides with the position of the two central galaxies 



the cluster probably had time to reach equilibrium after the 
merging event. 

The galaxy velocities have been transformed to the ref- 
erence frame of the cluster and in order to calculate the 
distances we have transformed the cluster velocities to the 
reference frame of the cosmic microwave background and 
assumed the ACDM cosmological model with parameters 
Qm = 0.3, = 0.7 and h = 0.7. The line-of-sight veloci- 
ties of 272 galaxies as a function of their projected distance 
from the cluster centre are plotted in the upper left panel 
of Fig. [1] We will refer to this kind of plot as a velocity dia- 
gram. We can see that although the main body of the cluster 
is well visible in the diagram, there is also a strong contam- 
ination by foreground and background galaxies which are 
unlikely members of the cluster. 

The third panel in the left column of Fig. [T] shows the 
velocity dispersion profile crios calculated from the sample 
with the standard unbiased dispersion estimator and n = 30 
galaxies per bin. The data p oints were assigned sampling er- 
rors of size (Tios/ •\/2(n — 1). We can see that the velocity dis- 
persion profile increases strongly with the projected distance 
R contrary to what is expected for a relaxed galaxy cluster. 
For such objects, real or simulated, the density profile is close 
to NFW (Navarro, Frenk & White 1997) while orbits do not 
depart significantly from isotropic and the line-of-sight ve- 
locity dispersion profile (7ios(i?) should be a decreasing func- 
tion of R (Lokas & Mamon 2001). We can try to reproduce 
this (Jios(-R) with the solution of the Jeans equation (Binney 
& Mamon 1982) 

where p{r) is the NFW density profile and T.m{R) its two- 
dimensional projection, assuming the anisotropy parameter 
(3 — 1 — ag{r)/ar{r) = and estimating the virial mass 
M„ and concentration parameter c (we assume c > 1). The 
results of the fitting procedure are shown in the second panel 
in the left column of Fig.[l]in terms of the 68.3, 95 and 99.73 
percent probability contours in the M„ — c parameter plane. 
We see that the preferred concentration is the lowest possible 
(c = 1) and the virial mass is as high as M„ = 2.3 x 10^^ 
Mq. The solution is plotted as a dashed line in the third 
panel of the left column of the Figure. 

The final modelling of galaxy cluster kinematics should 
take into account possible departures from isotropy of galac- 
tic orbits. Since the velocity dispersion data alone cannot 
break the degeneracy between the anisotropy and the shape 
(concentration) of the density profile it is useful to consider 
a higher order velocity moment, the line-of-sight kurtosis 
K-iosiR) = <I(^?)/o-L(i?) (Lokas & Mamon 2003; Lokas 
et al. 2006a). The kurtosis profile expressed in terms of 
k = [log (3if/2.68)]^'^^'' (where K is the standard kurtosis 
estimator) for the total sample of galaxies in A576 is shown 
in the lower left panel of Fig. [1] Note that contrary to the 
common belief that the higher moments of velocity should 
be more affected by interlopers (which is indeed often the 
case) the kurtosis profile is quite fiat and similar to profiles 
for regular clusters with very low contamination (see Fig. 5 
in Lokas et al. 2006a). 

Performing a joint fitting of velocity dispersion and kur- 
tosis profiles to the solutions of the Jeans equations (Lokas & 
Mamon 2003; Lokas et al. 2006a) assuming /3 = const we ar- 
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Figure 1. First row: velocity diagrams for the initial sample of galaxies (left) and after removal of interlopers with the method of Umax 
(middle) and the method using Mp/My^ (right). Rejected galaxies are marked with open circles, those included in the sample with 
filled circles. Second row: results of fitting ai^s assuming /3 = in the form of 68.3, 95.4 and 99.73 percent probability contours in the 
Mv — c parameter plane. The best-fitting values of the parameters are indicated with a dot. Last two rows: velocity dispersion and 
kurtosis parameter profiles for the corresponding samples of galaxies. Dashed lines show the best-fitting (Ti^g profiles obtained by fitting 
fTloa data alone assuming /3 = 0, solid lines show profiles from the joint fitting of ai^s and Kioa for arbitrary orbits. 



rive at an even higher mass estimate of Af„ — 4.0 x 10 Mq 
(again with c = 1) and anisotropy parameter (3 = —0.75. 
The parameters are summarized together with la error bars 
in Table[T]in the lines labelled as 'cut-off' referring to the ve- 
locity cut-off by which the galaxy sample was selected. The 
best-fitting solutions are plotted in the lower left panels of 
Fig.[T]as solid lines. The reason for the higher mass estimate 
in the case of fitting both velocity moments is the following: 
since kurtosis is not strongly overestimated it allows tan- 
gential orbits (see Fig. 10 in Sanchis, Lokas & Mamon 2004) 
which also fit the dispersion profile better, but for a higher 
mass. 



2.2 Removal of interlopers 

The method of interloper removal found to work best on 
simulated data (see Table 1 of Wojtak et al. 2006) was pro- 
posed by den Hartog & Katgert (1996). This approach was 
successful in removing on average 73 percent of unbound 
particles from simulated velocity diagrams obtained with the 
same initial cut-off in line-of-sight velocity, as applied here 
to A576: ±4000 km s~^ with respect to the cluster mean. 
The unbound particles remaining in the sample afterwards 
were those falling within the cluster main body in the pro- 
jected phase space. Those remaining interlopers do not bias 
the velocity dispersion significantly. 
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Table 1. Fitted parameters of A576 with Icr error bars. The first 
three rows show the results obtained from fitting only ctios for 
the samples obtained with simple velocity cut-off, the method 
of ^max and the method based on Mp/Myj' assuming isotropic 
orbits. The last three rows list the parameters obtained from joint 
fitting of (Tiog and /tjog for arbitrary orbits. 

method /3 M^IO'^'^Mq] c x^/N 

cut-off 23+5 l-O-o.o 9-4/'^ 

Mp/MvT 8-4li:6 10.21^^7° 9.2/5 

cut-off -0.75+°;55 40+2^ l-O^o.o 10.6/15 

0.1ll?:«t 9.01^:1 8.6li«gO 15.6/11 

Mp/MvT 0.07^?;^^ S-^llo lO-Ol^^'o" 16.4/11 



in 




r [Mpc] 

Figure 2. Mass profiles estimated from eq. Q in the iterative 
procedure of interloper removal based on rejection of galaxies with 
velocities larger than Umax given by eq. l[2]l. The highest line is 
for the first iteration, the lowest for the last. 



The method relies on calculating the maximum veloc- 
ity available to cluster members. The velocity is estimated 
assuming that a galaxy is either on a circular orbit with ve- 
locity Heir ~ \/ GM{r)/r or infalling into the cluster with 
velocity winf — v^Hcir so that 



max_R{uinf cos 9, v^r sin 9}, 



(2) 



where 9 is the angle between position vector of the galaxy 
with respect to the cluster centre and the line of sight. 

In order to calculate the maximum velocity we need an 
estimate of the mass profile. It turns out (see Wojtak et 
al. 2006) that the method works best if this is done with 
the mass estimator Mvt derived from the virial theorem 
(Heisler, Tremaine & BahcaU 1985) 



MvT(r = Rm^^) 



3ttN Eijvi -«)^ 
2G I]i< j 



(3) 



where A'^ is a number of galaxies with R < i?max, Vi is the 
velocity of the i-th galaxy and Rij is a projected distance 
between i-th and j-th galaxy. The mass profile can be then 
obtained from M{r) ~ AlyriRi < r < Ri+i), where Ri is 
the sequence of projected radii of galaxies in the increasing 
order up to 7?max = 3 Mpc in our case. From this mass profile 
we calculate the maximum velocity profile Q and remove 
galaxies with velocities exceeding Wmax. The procedure is 
repeated until no more interlopers are removed. 

The performance of the method in the case of A576 is 
illustrated in Fig. [2] showing the mass profiles in subsequent 
iterations with the highest line corresponding to the first 
iteration and the lowest to the last one. While in the first it- 
eration the mass within i?max = 3 Mpc is as high as 4 x 10^^ 
Mq, in the last iteration it decreases down to 1.5 x 10^^ ^Q- 
The resulting final sample of galaxies is shown in the upper 
middle panel of Fig. [1] as filled circles, while the rejected 
galaxies were marked as open circles. The corresponding ve- 
locity dispersion profile is plotted in the third middle panel 
of Fig. [1] (we restrict the analysis to seven inner data points 
so that all contributing galaxies are within the estimated 
virial radius) . It is now a decreasing function of the projected 
radius R, as expected for samples free of contamination. The 



best-fitting NFW parameters obtained from fitting aios are 
now My = 8.9 X 10^"' Mq and c = 9.2 with the confidence 
regions shown in the second middle panel of Fig. [1] 

Allowing for arbitrary constant anisotropy and adding 
kurtosis to the analysis (lower middle panel of Fig. [ij we 
arrive at very similar numbers, in particular the estimated 
anisotropy parameter is very close to zero (see Table [l|. 
The best-fitting profiles of velocity moments in this case 
are shown in the middle column of Fig. [l] with solid lines. 
Note that the la error bars estimated from joint fitting of 
both velocity moments, but keeping /3 as a free parameter, 
are larger than those found from fitting aios alone assuming 
/? = 0. 

Another method of interloper removal that we find use- 
ful in the case of A576 is the one based on the use of mass es- 
timators (see Table 1 of Wojtak et al. 2006). In this method, 
originally proposed by Perea et al. (1990), in addition to the 
virial mass estimator ^ we consider the projected mass 
estimator which for /3 = reads 



Mp = 



32 
TrGiV 



vfR^ 



(4) 



where again Vi and Ri are velocity and projected radius of 
the i-th galaxy and A'' is the number of galaxies in the sam- 
ple. As discussed by Wojtak et al. (2006) the presence of 
interlopers causes strong overestimation of both Mvt and 
Mp, but more so in the case of Mp, one can therefore con- 
struct a method of interloper removal based on the ratio 
Mp/MvT- Such a method is able to remove on average 65 
percent of unbound galaxies from velocity diagrams. 

The method uses the jackknife statistics in order to 
eliminate the interlopers. The data are arranged in {Ri,Vi} 
sequence, where i goes from 1 to n and we calculate n values 
of both mass estimators corresponding to n subsequences 
with one data point excluded. We identify as an interloper 
the galaxy for which the corresponding subsequence is the 
source of the most discrepant value of one of the estimators. 
In the next step, the same procedure is applied to a new data 
set with n~l particles. The disadvantage of this procedure 
is that it does not converge and we have to define some cri- 
terion for stopping the algorithm. Since in the case of sam- 
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Figure 3. Upper panel: the virial mass estimator Mvt (solid 
line) and the projected mass estimator Mp (dashed line) as a 
function of the number of removed galaxies rtrcm. Lower panel: 
the ratio of the two mass estimators Mp / Myx as a function of 
?irem. Vertical lines in both panels mark nrem = 44 when the 
procedure should be stopped. 

pies cleaned from interlopers we should have Mp /Mvt ^ 1 
and J^M/M — > these would be the most obvious condi- 
tions. Unfortunately, in the case of simulated clusters they 
are never exactly fulfilled and more conservative relations 
have to be adopted. 

However, applying the algorithm to the kinematic sam- 
ple of A576 we discovered that the condition Mp /Mvt — > 1 
is actually fulfilled after some interlopers are removed. The 
performance of the procedure is illustrated in Fig. O where 
we plot the values of the two estimators (upper panel) and 
their ratio (lower panel) as a function of the number of re- 
moved galaxies rircm. We find that at rircm = 44 the ratio 
Mp/MvT drops to unity and remains at this level even if 
more galaxies are rejected. This moment is a clear signature 
that the sample has been cleaned of interlopers. 

The final sample of galaxies from this method after re- 
moval of Tircm = 44 galaxies is shown in the upper right panel 
of Fig.[T] The galaxies accepted by the procedure are plotted 
with filled circles, while those removed with open circles. The 
corresponding velocity dispersion profile and the constraints 
on the parameters from the solution of the Jeans equation 
with /3 = are also shown in the right column of the Figure. 
The best-fitting parameters are now A'/„ — 8.4 x 10^* Mq 
and c = 10.2. Adding kurtosis to the analysis (lower right 



panel of Fig. [1} we again arrive at very similar numbers (see 
Table[l} . The best-fitting profiles of velocity moments in this 
case are shown in the right column of Fig. [T] with solid lines. 

The parameters estimated from fitting velocity disper- 
sion and both moments for the samples obtained by the two 
methods are summarized in Table [T] The present final sam- 
ple differs from the one obtained with the previous method 
only by one galaxy which is removed here and was retained 
in the previous case. This single galaxy is however respon- 
sible for the difference between the velocity dispersion and 
kurtosis profiles in the two cases (we keep 30 galaxies per bin 
so the binning of the data is different). The fitted param- 
eters of the NFW mass distribution however remain quite 
similar. 

Interestingly, no other method studied by Wojtak et al. 
(2006) and shown to work reasonably well for simulated clus- 
ters, is applicable to A576. In particular, the commonly used 
3(T method originally proposed by Yahil & Vidal (1977) does 
not remove any galaxy from the sample, whether applied 
to a single bin or in radial bins of different size. Similarly, 
the method based on fitting the isotropic solution of the 
Jeans equation and rejecting galaxies outside the 3o-ios(J?) 
lines successfully applied to less contaminated clusters by 
Lokas et al. (2006a) does not work here (it removes only 
one galaxy). The reason is that in the case of A576 we do 
not have single outliers but rather a uniform background of 
interlopers. This makes the dispersion profile for the con- 
taminated sample increase so strongly (see the left column 
of Fig.[T} that all galaxies fall within the fitted 3o"ios(i?) lines 
and the procedure does not start. 

2.3 Comparison with earlier work 

We studied the performance of a few methods of interloper 
removal, two of which work especially well: one based on esti- 
mating the maximum velocity available for cluster members 
and the other on the use of the ratio of two mass estima- 
tors. The methods work equally well producing reliable final 
galaxy samples that lead to similar conclusions about the 
mass distribution in the cluster. If we were to recommend 
one of these methods, we would choose the first since on av- 
erage it should remove more unbound galaxies and is always 
convergent. 

Other methods of selecting the cluster members have 
been applied to A576 before. For example, Mohr et al. (1996) 
used galaxies without emission lines as likely virialized mem- 
bers. However, the velocity dispersion profile for this sam- 
ple (see the upper panel of their fig. 8) is still rather fiat at 
larger distances from the cluster's centre leading to an over- 
estimated mass of 1.5 X 10^^ Mq (while the disp ersion of 
the 20 brightest members of this sample gives a value twice 
smaller) . This shows that selecting early-type galaxies is not 
sufficient to remove interlopers since some of these galaxies 
can still come from the cluster's neighbourhood. Rines et al. 
(2000) used the caustic method to study the infall patterns 
around A576 using a similar initial kinematic sample as the 
one used here. Their selection of members yields a similar 
sample as we obtain using our two methods of interloper 
removal. 

Since we believe the Umax method to be more reliable 
we adopt its results as final: M„ = 9.012 2 ^ 10^'' ^0 ^'^'^ 
c = 8.6lg*g'^. Our mass estimate is in excellent agreement 
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with the resuh of Rines et aL (2000). Within our virial ra- 
dius of 2.5 Mpc they find (see their fig. 7) the projected mass 
Mp = (10.1 ± 1.1) X 10^"* Mq which for our best-fitting con- 
centration corresponds to Af„ = (8.5 ± 1.0) x 10^* Mq. We 
also confirm their finding that the mass determined directly 
from traditional mass estimator Mvt strongly overestimates 
the true value. As shown in Fig. [2] the last iteration (lowest 
line) gives at 2.5 Mpc Mvt = 14.5 x 10^'' Mq which is 1.6 
times higher than our best estimate from velocity moments. 
The same is true for the projected mass estimator Mp which 
at the moment of stopping the interloper removal procedure 
is equal to Mvt (see Fig. |3]). In agreement with Biviano et 
al. (2006) we therefore find that these mass estimators over- 
estimate the cluster mass even if the interlopers are properly 
removed. 

Comparison with mass estimates derived from the con- 
dition of hydrostatic equilibrium of the X-ray emitting gas 
is more complicated because of a large number of parame- 
ters involved in such analyses. Using the result for Mx{< r) 
from Mohr et al. (1996, their equation 5.3) with T = 4 
keV (Kempner & David 2004) and r = 2.5 Mpc we get 
Mx{< 2.5) « 7 X 10^** Mq, which is consistent with our 
mass. (Restricting the analysis to r = 0.7 Mpc to avoid the 
extrapolation of X-ray data we get Mx{< 0.7) « 2 x 10^* 
Mq while our best estimate of mass within this radius is 
M = 3.4 X 10" Mq.) However, more detailed modelling of 
the X-ray data by Rines et al. (2000) showed that the mass 
profile inferred from the properties of X-ray gas is systemat- 
ically lower than the one derived from galaxy kinematics, as 
is the case for many clusters. Recently, Benatov et al. (2006) 
proposed that this discrepancy can be reduced if the galaxy 
orbits are radially anisotropic. Indeed, fixing /3 at the upper 
limit of the la confidence regions listed in the last two lines 
of Table[T]we find that the best-fitting M„ and c are pushed 
towards the lower limits of the ranges given in the Table. We 
therefore confirm that significantly lower mass estimates are 
preferred for radial orbits. 



3 OTHER CLUSTERS 

We have shown that for strongly contaminated kinematic 
samples, as in the case of A576, the proper removal of inter- 
lopers is essential to avoid bias in the estimated parameters 
of the mass distribution. One may ask to what extent A576 
is special and whether we should worry about such cases or 
perhaps just discard them from our cluster samples. As dis- 
cussed in Lokas et al. (2006a) there are now few tens nearby 
{z < 0.1) clusters with significant number of galaxy red- 
shifts available from NED. About half of these clusters have 
irregular X-ray images so are far from equilibrium and there- 
fore are not tractable by Jeans modelling. In addition, some 
of the clusters, is spite of regular X-ray images, have close 
neighbours or strong substructure which makes the analysis 
difficult. A576 is a special case among these clusters in a 
sense that it has a strong contamination of interlopers but 
otherwise is quite regular. 

A576 is not the only cluster which can be retained for 
dynamical analysis after careful interloper removal. To il- 
lustrate further the importance of interloper removal and 
show broader impact of the schemes studied in Wojtak et 
al. (2006) we studied 15 more nearby (z < 0.1) and possibly 



Table 2. Properties of the cluster sample composed of 15 clusters. 
The columns list the cluster name, redshift, total number of galax- 
ies A^tot in the initial sample and numbers of interlopers identified 
with three methods of interloper removal. The boldfaced numbers 
indicate interlopers excluded for the final modelling. The last col- 
umn gives the main sources of the redshift data: 1 - Durret et al. 

(1998) , 2 - Katgert et al. (1998), 3 - Smith et al. (2004), 4 - 2dF 
Galaxy Redshift Survey (CoUess et al. 2001), 5 - SDSS (Abazajian 
et al. 2005), 6 - SDSS (Abazajian et al. 2004), 7 - Chincarini & 
Rood (1971), 8 - Lucey et al. (1997), 9 - Scodeggio et al. (1995), 
10 - Ferrari et al. (1998), 11 - Pimbblet et al. (2006), 12 - Hill & 
Oegerle (1993), 13 - Barmby & Huchra (1998), 14 - Wegner et al. 

(1999) . 
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relaxed (in terms of regularity of their X-ray maps and ve- 
locity diagrams) clusters with at least 120 galaxies available 
within projected radius of 2 Mpc. We apply to all of them 
different procedures of interloper removal and compare cor- 
responding samples of cluster galaxies in terms of numbers 
of rejected galaxies and results of fitting the isotropic solu- 
tion of the Jeans equation. For clarity we focus on three most 
representative methods described in detail in the previous 
sections: 3a{R) rejection, Wmax and Mp/MvT (see Wojtak 
et al. 2006), since other approaches are very similar. Like 
for A576, we use NED to select positions and redshifts of 
galaxies. Projected distances and the line-of-sight velocities 
of galaxies were obtained following the same steps as de- 
scribed in section 2 for A576. Table [2] lists the redshifts, the 
total numbers of galaxies used and the results of interloper 
removal for the 15 selected clusters. The Table also gives the 
references to the most important sources of redshift data. 

We have grouped the clusters in Table [5] according to 
their behaviour in terms of interloper removal. The results 
of interloper removal for the first five clusters in Table [2] are 
illustrated in Fig. ldal In the first three columns of Fig. Idal we 
put velocity diagrams with filled and open circles indicating 
cluster galaxies and identified interlopers. Each row corre- 
sponds to a different cluster as marked in the Figure. The 
columns from the left to the right correspond respectively 
to the results of 3a{R), Wmax and Mp/Mvt procedures. The 
last column shows results of fitting isotropic solution of the 
Jeans equation, ^ with /3 = 0, to the velocity dispersion 
profiles in the form of 68.3, 95.4 and 99.73 percent prob- 
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Figure 4a. Results of interloper removal for the first group of clusters in Table [2] Three first columns from the left to the right show 
results of the three methods of interloper removal: 3cr{R), i)max and Mp/MvT, respectively. Filled circles indicate cluster members and 
open circles the galaxies identified as interlopers. The right column shows results of fitting ctios to the solution of the Jeans equation JTJ 
with /3 = in terms of 68.3, 95.4 and 99.73 percent probability contours in the My — c parameter plane. Dotted, solid and dashed lines 
correspond to samples of galaxies obtained with 3cr(/?), Umax and Mp/Myx procedures of interloper removal. In the cases where only 
two types of contours are shown in the right column, two methods gave identical results. 



ability contours in Mv — c plane. Dotted, solid and dashed 
lines correspond to samples of cluster galaxies obtained with 
3a{R), fmax and Mp/Mvt methods. Similar plots for the 
two other groups of clusters listed in Table [5] are shown in 
Fig.liblandlicl 

Since each cluster is different it is difficult to divide 
them into distinct classes from the point of view of their 
behaviour under interloper removal schemes. However, com- 
paring the numbers of interlopers determined with the three 
methods (Table [2j we find that for 5 out of 15 clusters (the 
first group in Table [2]) all results of interloper removal are 
nearly the same. A good example of this class of clusters is 
the A85 shown in the first line of plots in Fig. I4al We notice 
that in this case only the Mp /Mvt procedure seems to give 



slightly less consistent output. This happens due to the fact 
that the number of interlopers here is relatively small and 
it is hard to find precisely a characteristic knee-like point 
in the Mp/Mvt profile (see Fig. |3]for A576) at which the 
algorithm of interloper removal should be stopped (Wojtak 
et al. 2006) . However, the differences between the final sam- 
ples of cluster galaxies are not significant and, consequently, 
probability contours in the M„ — c plane determined from 
the data sets overlap each other. 

Five clusters of the second group in Table [2] illustrate 
the case when only i;max and Mp/Mvr schemes deter- 
mine samples of cluster galaxies consistent with each other, 
whereas the 3cr(ii) approach fails giving rather discrepant 
results. This situation manifests itself clearly in terms of 
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Figure 4b. The same as Fig. l4al for the second group of clusters in Table [2] 



the fitting results. Probability contours associated with Umax 
and Mp/MvT schemes coincide almost exactly, whereas the 
ones corresponding to the "ioi^R) procedure are shifted to- 
wards higher mass and lower concentration. Note that from 
the point of view of interloper removal A576 is very similar 
to the clusters of this group. 

Finally, for the last 5 clusters listed in Table [5] all meth- 
ods of interloper removal identify different numbers of in- 
terlopers. This implies that probability contours associated 
with the selected samples of cluster galaxies typically do not 
coincide as well as in the case of clusters of the first section 
of Table [2] Among these clusters A2147 and A4038 have 
distinct neighbouring substructures which overlap partially 
with the main body of the cluster. This causes the failure of 
both Sa{R) and Vmax procedures. Nevertheless, the results 
obtained with the Mp/Mvt scheme are likely reliable for 
these two clusters since the moment of convergence of the 
algorithm is clearly seen in the Mp/Mvt profile. 



4 DISCUSSION 

The number of interlopers in a sample for a given cluster 
depends on its particular neighbourhood but also on the 
method of galaxy selection by the observer (whether some 
photometric criteria were used). The fraction of interlopers 
can therefore be very different and is not so important. More 
important is their velocity distribution. If there are only a 
few interlopers but with highly discrepant velocities they can 
also strongly bias the inferred parameters of mass distribu- 
tion. However, such interlopers are typically easily removed 
by simple methods like the rejection of the 3cr outliers. In 
A576 the interlopers are more uniformly distributed in ve- 
locity which requires the use of more sophisticated meth- 
ods like the ones presented here (still, the methods are sim- 
pler in application than the caustic method of Rines et al. 
2000). Since regular clusters are scarce and the problem of 
contamination by interlopers is common among them we 
think even those most contaminated should be studied. We 
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Figure 4c. The same as Fig. I4al for tlie last group of clusters in Table[2] 



demonstrated that in many such cases proper removal of in- 
terlopers can lead to reliable determination of galaxy mem- 
bership. 

The general picture that emerges from our study of 16 
clusters is that even in the case of significant contamination 
different interloper removal schemes can yield similar results 
increasing our confidence in the final galaxy samples. In the 
case when two or more schemes (usually the methods based 
on Umax and Mp /Mvt) give similar results we recommend 
the use of Umax since it converges and is independent of any 
assumptions. There are cases however with strong substruc- 
ture visible in the velocity diagrams when only the Mp /Mvt 
method seems to give results which clearly separate the sub- 
clusters from the main body of the cluster of interest. We 
recommend this method for such cases. 

Following these conclusions we have selected final 
galaxy samples for our 15 clusters presented in Fig. 4a-c 
of the previous section: we used the Mp/Mvr method for 
A2147 and A4038 and the «max method for all the remain- 



ing clusters. The resulting profiles of the velocity moments 
are shown in Fig. [S] and |S] As for A576, we fitted the data 
for the velocity dispersion and kurtosis to the solutions of 
the Jeans equations with /3 = const estimating the virial 
mass, concentration and anisotropy parameter. The best- 
fitting parameters together with la error bars are given in 
Table [S] We have also listed there for comparison the pa- 
rameters obtained previously for the same galaxy samples 
from fitting velocity dispersion alone assuming /3 = 0. The 
fitted profiles are shown in Fig. [5] and |6] with solid lines (ve- 
locity dispersion profiles obtained for /3 = are plotted in 
Fig. [5] with dashed lines which in most cases overlap almost 
exactly with the solid lines). 

These new results allow us to extend the mass- 
concentration relation given in Lokas et al. (2006a) to the to- 
tal of 22 clusters. The relation is shown in Fig. [7]with circles 
plotting the parameters for the 6 clusters studied in Lokas 
et al. and squares for the 16 new clusters studied here. The 
dashed line shows the approximation of the c{Mv) relation 
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Figure 5. Profiles of the velocity dispersion for 15 clusters listed in Table [2] In all cases the Dmax method of interloper removal was 
applied except for A2147 and A4038 where the method based on Mp/MvT was used. The dashed lines plot the best-fitting velocity 
dispersion profiles obtained for /3 = 0, the solid lines plot the solutions obtained when fitting both dispersion and kurtosis. The best-fitting 
parameters are listed in Table [3] 
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Figure 6. Profiles of the kurtosis variable k = [log (3i<'/2.68)]^'^^'^ for the 15 clusters listed in Table[2] The galaxy samples used are the 
same as for Fig. [5] The solid lines plot the solutions of the Jeans equations obtained when fitting both dispersion and kurtosis. 



calculated from the toy model proposed by Bullock et al. 
(2001) which reproduces well the properties of a large sam- 
ple of haloes found in their A''-body simulations. The fit of 
a linear relation of the form log c = a log[M„/(10^'*MQ)] + h 
to the data gives best-fitting parameters significantly dif- 
ferent from those in Lokas et al.: a = —0.17 ± 0.65 and 



h = 1.14 ± 0.62 (at 68 percent confidence level). This best- 
fitting relation is shown in Fig. [7] as a solid line: the slope of 
the fitted relation agrees well with the one from A^-body sim- 
ulations, but our concentrations are shifted towards higher 
values. The results agree within errors with other determina- 
tions of the mass-concentration relation (e.g. Pointecouteau, 
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Table 3. Fitted parameters of 15 clusters from Table [2] with Icr 
error bars. The parameters were obtained from fitting ctios with 
[3 = and joint fitting of trios and ki^s for arbitrary constant 
/3. The galaxy samples for A2147 and A4038 were obtained with 
the Mp / Myt method, for the rest of the clusters with the rimax 
method. 
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Figure 7. The mass-concentration relation. The points show the 
best-fitting My-c pairs for 22 clusters with la error bars. The 
circles correspond to the 6 clusters studied in Lokas et al. {2006a), 
the squares to the 16 new clusters studied here. The dashed line 
plots the prediction from the A'-body simulations by Bullock et 
al. (2001) while the solid line shows the best fit to the data points. 



Arnaud & Pratt 2005; Vikhlinin et al. 2006; Mandelbaum 
et al. 2006; Rines & Diaferio 2006; Schmidt & Allen 2006; 
Buote et al. 2006). Interestingly, our higher typical concen- 
tration values agree with other recent determinations for 
individual clusters (e.g. Broadhurst et al. 2005) . 

Note that the profiles of the velocity moments in Fig. [5] 
and [6] determined from the cleaned galaxy samples are very 
regular, with expected behaviour: the dispersion profile de- 
creases slowly at larger distances while the kurtosis profile 
is rather fiat. The profiles are well behaved even in the 
case of A119, A2634 and A2147 in spite of the fact the 
surface brightness distributions of the X-ray gas are least 
regular in these clusters. From this point of view there is 
no difference between these originally contaminated clusters 
and those presumably more regular ones studied in Lokas 
et al. (2006a). The difference between these two groups 
manifests itself however in the fitted parameters. While in 
the previous study 4 out of 6 clusters had the best-fitting 
anisotropy parameter |/3| < 0.2, for the present sample this 
is only the case for 5 clusters out of 16. Two of the clus- 
ters (A1651 and A2670) have -2 < /? < -1.6 correspond- 
ing to ae/ffr = 1.6 — 1.7, indicating significant tangential 
anisotropy, while three clusters (A85, A3571 and A4059) 
have (3 « 0.7 corresponding to drjog = 1.8, indicating sig- 
nificant radial anisotropy. This may point towards a more 
perturbed state of the present sample. Still, all clusters have 
fi consistent with isotropy at Icr level, except for A85. 

A85 seems to be the most suspicious case in the present 
sample. While the removal of interlopers would suggest no 
problems in the analysis (all three methods yield very simi- 
lar results in terms of the number of rejected interlopers as 
well as the fitted mass and concentration, see Table [2] and 
Fig. I4ap . while fitting both velocity moments we arrive at 
rather different parameters, lower mass and concentration 
and more radial orbits. The profiles of the velocity moments 
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are also not typical: the first two dispersion data points are 
much higher than the rest and the kurtosis points lower. 
This may point at departures from equilibrium which are 
indeed confirmed by Kempner, Sarazin & Richer (2002) and 
Durret, Lima Neto & Forman (2005) who studied the X- 
ray data from Chandra and XMM-Newton and discovered 
a subcluster about to merge with A85 as well as evidence 
for past merging events. Note, however, that the mass esti- 
mated from X-ray data by Durret et al. within r < 1.3 Mpc, 
Mx ~ 4.3 X 10^*Mq, is in very good agreement with our 
mass within this radius M — 4.5 x 10^'' Mq. 

Reiprich & Bohringer (2002) used the ROSAT X-ray 
data to estimate the masses of all clusters studied here ex- 
cept for A2593 and A2670. In all cases our estimates agree 
with their findings within errors, only for A426, A3667 and 
A4038 our masses are about a factor of two higher. Our 
mass estimates also agree within errors with the ones found 
from the virial theorem by Girardi et al. (1998) in the case of 
A85, A119, A401, A426, A576, A2593, A2634, A2670, A3571 
and A3667. For A2142 our mass estimate is more than twice 
lower than the one of Girardi et al. The result in this case 
may be strongly dependent on the sample selection because 
the velocity diagram of the cluster is rather irregular. 

For A3921 our estimate of the virial mass is almost four 
times higher than the one of Girardi et al. which is prob- 
ably due to their rather small sample of only 29 member 
galaxies while our sample contains 197 members. Referring 
to X-ray studies with XMM-Newton (Belsole et al. 2005) 
we find that there are signatures of recent merger in the 
form of a bar-like structure of enhanced temperature. Bel- 
sole et al. (2005) estimate the mass to be Mx = 2 x IO^Mq 
within r < 0.8 Mpc while within this distance we get 
M = 4.2 X 10^* Mq, this time a value a factor of two higher. 
Note however, that our estimate is in very good agreement 
with the older ROSAT data (Reiprich & Bohringer 2002) 
which give M„ = 11 x IO^'^A/q. Going back to the velocity 
diagram for A3921 in Fig. I4bl we see that there are some 
doubtful galaxies (e.g. the chain at 2000 km s~^) which are 
not rejected by our algorithms but which may not be clus- 
ter members and may in fact overestimate the true cluster 
mass. 

A2147 is a member of the Hercules supercluster which 
obviously makes the assignment of galaxies to each clus- 
ter more difficult (Tarenghi et al. 1980; Barmby & Huchra 
1998). Together with A4038 it is the only cluster for which 
neighbours are well visible and we believe the Mp/Myr 
method of interloper removal is more reliable than Umax- In 
the case of A2147 our mass estimate is about twice lower 
than the one based on velocity dispersion by Barmby & 
Huchra (1998). The difference is obviously due to the dif- 
ferent selection of galaxies, for example they did not remove 
the group at —1300 km s~^ (see Fig. l4c)l . In the case of 
A4038 we agree with the rough estimate of the cluster mass 
based on galaxies in the core by Green, Godwin & Peach 
(1990). 

We have demonstrated that after careful removal of in- 
terlopers a significant number of galaxy clusters can be made 
eligible for Jeans analysis. Using the example of A576 we 
have shown in detail how the different most reliable meth- 
ods of cleaning the kinematic samples of interlopers work 
on real data. We illustrated the complexity of the problem 
with examples of other clusters displaying different degrees 



of contamination and different distributions of interlopers. 
The results of the final dynamical analysis of these clusters 
are not significantly different from those obtained for more 
regular clusters increasing our confidence in the applicability 
of Jeans formalism. 
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